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We propose to study the transport through tunneling-coupled double quantum dots (DQDs) 
connected in series to leads, using the fmite-(7 slave-boson mean field approach developed initially 
by Kotliar and Ruckenstein [Phys. Rev. Lett. 57, 1362 (1986)]. This approach treats the dot-lead 
coupling and the inter-dot tunnelling t nonperturbatively at arbitrary Coulomb correlation U, thus 
allows the anti- ferromagnetic exchange coupling parameter J = 4t 2 /U to appear naturally. We find 
that, with increasing the inter-dot hopping, the DQDs manifest three distinct physical scenarios: 
the Kondo singlet state of each dot with its adjacent lead, the spin singlet state consisting of local 
spins on each dot and the doubly occupied bonding orbital of the coupled dots. The three states 
exhibit remarkably distinct behavior in transmission spectrum, linear and differential conductance 
and their magnetic-field dependence. Theoretical predictions agree with numerical renormalization 
group and Lanczos calculations, and some of them have been observed in recent experiments. 
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The prediction of the Kondo effect in a quantum dot 
has been verified by observations, m Very recently, with 
the experimental breakthroughs, much theoretical 
interest has been payed to transport through tunneling- 
coupled double quantum dots (DQDs) to investigate 
the interplay between the Kondo effect and the anti- 
ferromagnetic correlation. Aono et al |4|J^] and Georges 
et al |6| employed an infinite-t/ slave-boson mean-field 
(SBMF) theory to the problem. Aguado and Langreth 
generalized this method to an out-of-equilibrium sit- 
uation. However, the predicted abrupt jump in conduc- 
tance at certain condition, was proved to be an artifact of 
the model. || Furthermore, the necessity of the artificial 
introduction of a magnetic coupling parameter indicates 
the weakness of this kind of infinite-?/ SBMF theory in 
treating the kinetic exchange process in DQDs. J§| 

In this letter we propose to study the transport 
through DQDs using the finite-// SBMF approach, ini- 
tially developed by Kotliar and Ruckenstein for other 
system, |)]-[llj which enforces the correspondence be- 
tween the impurity fermions and the auxiliary bosons 
to a mean-field level to release the U — oo restriction. 
In DQDs this approach allows one to treat both the dot- 
lead coupling and the inter-dot tunnelling t nonperturba- 
tively at arbitrary strength of Coulomb correlation U and 
thus an anti-ferromagnetic exchange coupling parameter 
J = 4t 2 /U to appear naturally. The derived conduc- 
tance exhibits smooth change with changing tunneling 
strength and magnetic field rather than the abrupt tran- 
sition showing up in the previous SBMF theory. The 
agreement of the present results with those of numeri- 
cal renormalization group (NRG) and numerical Lanczos 
calculation justifies the suitability of the proposed finite- 
ly SBMF approach to DQDs. 

Transport through tunneling-coupled double quantum 
dots L and R (assuming a single level e n in the 77th dot, 



7] = L,R) connected separately to leads L and R (hav- 
ing energy e v k for Bloch state k, r\ = L, R) can be mod- 
eled by a two-impurity Anderson Hamiltonian with intra- 
dot Coulomb repulsion U v and dot-lead coupling V v , and 
inter-dot hopping t. In the presence of a magnetic field 
B the electron energies are split into e r)cr = e, ; + ah-Q and 
f-rika = f-r)k + c/ib for up (<t = 1) and down (<7 = — 1) 
spins, where 2/ib = 5MB-B (<? is the Lande factor and fj,B 
the Bohr magneton) . The physical problem of DQDs de- 
scribed by this Hamiltonian can be dealt with by means 
of the finite-?/ slave-boson technique of Kotliar and Ruck- 
enstein |§-(ll]] using the following effective Hamiltonian 
in an enlarged Hilbert space which contains, in addition 
to the electron (fcrmion) annihilation (creation) opera- 
tors c v ka ( c lfka) m t ne leads and c va (cj^) in the dots, a 
set of four auxiliary boson operators e v (e^), p va (pj^.) 
and drj (djj) for the 77th dot, which act, respectively, as 
projectors onto the empty, singly occupied (with spin up 
and down) and doubly occupied electron states: 
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Here, the operator z va is given as 
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The constraints, i.e. the completeness relation 
'YjoPncrVri*' + e v e v + dld^ — 1 and the condition for the 
correspondence between fermions and bosons c\ ltJ c 7la — 
Pno-Pw* + dtd v , which should be imposed for each dot, 
have been incorporated with Lagrange multipliers A^ 1 "* 

(2) 

and X^J in the effective Hamiltonian. The form (2) of 
z v<J operators is chosen to reproduce the correct U — > 
limit in the mean field approximation. 

From the effective Hamiltonian (1) one can derive 
equations of motion of the slave-boson operators. Then 
we use the mean-field approximation in these equations 
and in the constraints and the effective Hamiltonian, in 
which all the boson operators are replaced by their ex- 
pectation value. It is believed that at low temperatures, 
the mean field approximation is correct for describing 
the Kondo-type transport through QD in the Kondo 
regime. The Fermi operators in these equations are bi- 
linear and their expectation value can be expressed, with 
the help of the Langreth analytical continuation rules for 
the close time-path Green's function [[l3|, by the distri- 
bution Green's functions of both dots 
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(j] =/= rf), in which /^(w) is the Fermi distribution function 
of 77th lead, 

A a (uj) = [(U - Iho + iT La )(uJ - £Rcr + lT Ra) - 

X [(u> - l La - iTLa){u - £fl(T - iTro) - t%\, (4) 
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r r) = ^^2k l^>;| 2 ^( a -' ~ e vk) being the coupling constant 
between the 77th QD and the lead. These 8 equations, 
together with the statistical average of the 6 constrains, 
forms the closed set of 14 equations (ry = L.R and 
<7 = 1,-1) for total fourteen unknowns (e^, p V a, d^, 
and A^ 1 -*, The current through DQD is given by 

1 = f /-» r(w)[/L(w) - /hH], in which 
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is the transmission probability. 

This letter focuses on the situation of two identical 
QDs: e L = e R = e d , V L =Vr = V and T L = T R = T. In 
the linear limit the conductance can be written as 
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In the case of t <C U, the anti-ferromagnetic ex- 
change coupling can be described by J = 4t 2 /U. 



The Kondo temperature of the QD is given by Tk = 
C/ % /^exp(-7r//3)/2 7 r [(3 = -2UT / e d (U + e d )], where the 
energy level of the QD, e d , is measured from the Fermi 
level of adjacent lead. Hereafter in the text and in the 
figures we will always use T as the energy units. 

Fig. 1 shows the linear conductance G and the total 
occupied number N for t = 1,2 and 8 at U — 14, as 
functions of e d in the absence of magnetic field. Depend- 
ing on the hopping strength t, three different physical 
scenarios can be seen. (1) In the case of t < 1, when 
inter-dot coupling is weaker than the dot-lead coupling 
( J/7k 0.57) G has a single peak in the N ~ 2 regime 
(—0.75 < e d /U < —0.25), indicating that Kondo sin- 
glet states consisting of each individual dot and adjacent 
lead dominate. The inter-dot tunneling serves mainly as 
the conducting bridge across the two dots. The peak 
conductance increases with increasing t and reaches the 
unitary value 2e 2 /h at t = 1. (2) In the case of t = 2, 
when inter-dot coupling becomes stronger than the dot- 
lead coupling (J/Tk = 2.3), G descends slightly in the 
N ~ 2 range, exhibiting a shallow valley. The increas- 
ing inter-dot tunneling together with the strong intra-dot 
Coulomb repulsion favors the formation of a spin singlet 
state consisting of local spins on each dot. (3) In the 
case of t — 8, when the inter-dot correlation is strong 
enough to compensate the intra-dot Coulomb repulsion 
(t ^ U/2), the doubly occupied bonding (even) orbital of 
the coupled dots are favorite. G exhibits a double-peak 
structure with a deep valley almost down to zero, and 
the Kondo feature disappears completely. 

Fig. 2 shows the inter-dot-hopping-dependent conduc- 
tance. One can see the continuous appearance of the 
three different scenarios when the inter-dot hopping 
changes. G increases with increasing t from zero, reaches 
the maximum at around t ~ 1 and then decreases. The 
single peak structure of G at U — 8 turns into a pattern 
having a shoulder for larger [/, and finally the shoulder 
becomes a monotonously growing curve at U — > 00 limit. 
These results are in fairly good agreement with those of 
numerical Lanczos calculation [Q and NRG calculation 
H for finite Coulomb interaction U. 

In the present model the energies e e ( ) CT = e da =F t a can 
be considered as the energy levels of the even and odd 
states (in the case without the magnetic field they are 
independent of the spin and we will omit the a index). 
The energy separation between the odd and even states 
e — e e = 2t equals 2t for U = and is largely suppressed 
by the intra-dot Coulomb repulsion. On the other hand, 
in the Kondo-dominanted regime (J/Tk ^$ 1, N ~ 2) the 
resonant energy level e d locates very near to the chemical 
potential of the lead, such that the ratio P = e 2 d /(t/t) 2 is 
much smaller than 1 [Fig. 2(b)], and the conductance (6) 
is simplified to G — ^f- n^rp 1 which reaches its maxi- 
mum value at t = 1. 

The physical origins of the transport properties can be 
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understood in terms of the transmission spectra T(u>), 
which are depicted in Figs. 3(a)-(c) for = —2 and 
U = 8, 14, 20 in the absence of voltage bias. T(ui) reaches 
the maximum 1 at lu = e e and e Q . In the case of t = 1 
where the intra-dot Kondo correlation is predominant 
over the inter-dot coupling {J/T K = 0.78, 0.57, 0.43 for 
U = 8,14 and 20, respectively), e e and e are closed 
together that T(lu) exhibits only a single peak. In the 
case of t = 2 where the inter-dot coupling grows over 
the intra-dot Kondo correlation (J/Tk = 3.1, 2.3, 1.7 for 
U = 8, 14 and 20, respectively) two peaks in T(u>) show 
up [Fig. 3(b)]. An example of strong inter-dot coupling 
(t = 8) is plot in Fig. 3(c), where two peaks in T(u>) 
are so greatly separated that the valley bottom is nearly 
opaque. For = —2 the zero-frequency line is in be- 
tween e e and e D but always closer to the even peak. 

A more significant and experimentally measured quan- 
tity is the differential conductance dl/dV, which have 
been calculated at finite source-drain biases by assum- 
ing a symmetric voltage drop, [1l = —Hr = eV/2. The 
calculated results of dl/dV are plotted in Fig. 3(d)- (f) 
as a function of eV/Tn- The differential conductance 
behavior changes continuously with increasing i, but ex- 
hibits remarkably different characteristics for the above 
three typical scenarios. The predicted double peak of the 
Kondo resonance in the dl/dV-V curves has been ob- 
served in a recent measurement (in Ref. S). The sharp 
dips with negative differential conductivity showing up 
in the theoretical curves, may be smeared out by the de- 
phasing processes at finite bias voltages and other com- 
plexity in a real experiment. We look forward to further 
subtle experiment to test it. 

Fig. 4 shows the linear conductance in the presence of 
external magnetic fields (MC), as a function of Zeeman 
energy h% (a) and as a function of e^/U (b). Again, 
t = 1 and t = 2 represent two distinct physical situations. 
This behavior of the MC can be understood from the 
magneto-transmission spectra as shown in Figs. 5(a)-(c). 
In the case of t = 1, for which the T(u>) exhibits a single 
broad maximum around lu = at = 0, the trans- 
mission probability, and thus the linear conductance G, 
decline with increasing h-Q , indicating a negative MC and 
a weakening of the Kondo effect [Fig. 5(a)]. The stronger 
the intra-dot Coulomb repulsion, the faster the conduc- 
tance G decreases with increasing h-Q- When t > 2, both 
the even peak and odd peak in the transmission spec- 
trum T(lo) in the absence of the magnetic field may split 
into two peaks corresponding to up and down spins for 
strong magnetic fields, giving rise to a four-peak struc- 
ture [Fig. 5(c) for h-Q = 0.4]. Note that in the presence 
of intra-dot Coulomb repulsion the splitting between the 
up and down peaks is smaller than the Zeeman energy 
2/ib (similar to the situation of the single QD [0,0) and 
always wider in the even channel than in the odd channel. 

In the case of t = 2, where the even and odd peaks 
are moderately separated and the valley in T(oj) locates 



not far from lu = for = 0, the up-spin peak of the 
even channel may emerge with the down-spin peak of the 
odd channel at suitable magnetic fields, giving rise to a 
triple-peak pattern. As a consequence, T(lu) at lu = 
increases with increasing magnetic field up to a certain 
strength, yielding a positive MC and an enhancement of 
the Kondo effect. Of course, further splitting of the up 
and down spin peaks by the continuing increase of the 
magnetic field will lead to rapid drop of the transmission 
probability (and the conductance G) after reaching the 
maximum. As can be seen from Fig. 4(a) the position 
of the conductance maximum, is strongly dependent on 
the Coulomb interaction U. It should be noted that, this 
rapid drop, is apparently smooth rather than abrupt. A 
similar behavior of MC has been found by Aono et. al. 
Q in their U —> oo model, but their prediction that MC 
suddenly becomes zero at a certain magnetic field, seems 
to be an artifact of the model. 

Finally, we suggest that the magneto-transmission 
probability can be observed in experiment by measur- 
ing the dl/dV-V characteristic in the presence of mag- 
netic fields. Figs. 5(d)-(f) demonstrate that dl/dV tran- 
sits from single zero-bias peak to multiple-peak pattern 
depending on the inter-dot hopping t. 

In conclusion we have studied the transport through 
DQDs using the finite-?/ SBMF approach. Due to com- 
petition between the Kondo and anti-ferromagnetic cor- 
relations, the DQDs manifest three distinct physical sce- 
narios as the inter-dot hopping t increases from zero: 
the Kondo singlet state of each dot with its adjacent 
lead, the spin singlet state consisting of local spins on 
each dot, and the doubly occupied bonding orbital of 
the coupled dots. Although transition between them is 
continuous, the three states exhibit remarkably distinct 
behavior in transmission spectrum, linear and differential 
conductance and their magnetic-field dependence. These 
features have been partly observed recently, and further 
experiments are needed for a full test of the theory. 
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Figure Captions 

Fig.l Linear conductance G (thick lines) and total 
charge number N (thin lines) vs. the QD energy level 
€d at U = 14 for inter-dot hopping t = 1 (solid line), 2 
(dashed line), and 8 (dotted line) repectively. 

Fig. 2 (a) Inter-dot hopping-depcndent linear conduc- 
tance through DQDs with e<j = —2 at different Coulomb 
interaction U = 8, 14, 20, and oo. (b) The parameter 
P vs. the inter-dot hopping t for U = 14 and different 
energy levels €d = 0, —1, —2, —4, —5 and —6. 

Fig. 3 Transmission spectrum [(a)-(c)] and dl/dV-V 
curves [(d)-(f)] for = —2 and different values of the 
hopping t = 1, t = 2, and t = 8. 

Fig. 4 (a) The magneto-conductance vs. magnetic field 
Jib for Ed — —2. (b) The magneto-conductance vs. the 
QD energy level for hs = 0, 0.1, and 0.2. 

Fig. 5 Magneto-transmission spectrum [(a)-(c)] and the 
dl/dV-V curves [(d)-(f)] for U = 8, e d = -2 at different 
magnetic fields. 
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